Bose-Einstein condensation of stationary-light polaritons 
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We propose and analyze a mechanism for Bose-Einstein condensation of stationary dark-state 
polaritons. Dark-state polaritons (DSPs) are formed in the interaction of light with laser-driven 3- 
level A- type atoms and are the basis of phenomena such as electromagnetically induced transparency 
(EIT), ultra-slow and stored light. They have long intrinsic lifetimes and in a stationary set-up with 
two counterpropagating control fields of equal intensity have a 3D quadratic dispersion profile with 
variable effective mass. Since DSPs are bosons they can undergo a Bose-Einstein condensation at 
a critical temperature which can be many orders of magnitude larger than that of atoms. We show 
that thermalization of polaritons can occur via elastic collisions mediated by a resonantly enhanced 
optical Kerr nonlinearity on a time scale short compared to the decay time. Finally condensation 
can be observed by turning stationary into propagating polaritons and monitoring the emitted light. 
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Since the critical temperature of Bose-Einstein con- 
densation (BEG) pj, l2| is inversely proportional to the 
mass of the particles already early on in the history of 
the field quasi-particles where considered as candidates 
for condensation at high temperatures [3] . Very recently 
this subject regained a lot of attention triggered by exper- 
imental breakthroughs in microcavity exciton-polaritons 
y,|5|,|6[ and thin- film magnons [7]. Yet major drawbacks 
of these systems are the limitation to quasi-condensation 
in two dimensions, and in the case of excitons, the ne- 
cessity of a cavity which makes a distinction from lasing 
difficult [^. We here propose a mechanism for true con- 
densation in three spatial dimensions employing a differ- 
ent kind of polaritonic quasi-particles, called dark-state 
polaritons (DSP) [9]. DSPs |10j] emerge in the Raman 
interaction of lasers with 3-level quantum systems and 
are the basis of ultra-slow [11], stopped [l2|, ll3[, and sta- 
tionary light [14]. They have a considerably longer life- 
time than exciton polaritons or magnons, provide a 3D 
quadratic dispersion with variable mass, and have very 
high condensation temperatures. DSPs can easily be cre- 
ated and thermalization can be achieved on a time scale 
much shorter than their lifetime. Finally condensation 
can easily be observed by transforming stationary DSPs 
into light pulses. 

If an optically thick ensemble of three-level quantum 
systems is irradiated by a strong coherent coupling laser 
it can become transparent for a probe field within a cer- 
tain frequency range close to the two-photon Raman reso- 
nance. Associated with this phenomenon known as elec- 
tromagnetically induced transparency (EIT) [l5|, [l6| is 
the formation of polaritonic eigensolutions, called dark- 
state polaritons (DSP) [10|, which are superpositions of 
the probe-field amplitude and a collective Raman excita- 
tion. The mixing angle that determines the admixtures 
of field and matter degrees of freedom can be changed by 
varying the strength of the control laser. At the trans- 
parency frequency the dispersion of dark-state polaritons 
is linear with a slope determined by the mixing angle. 



This situation changes if two counter-propagating control 
laser with comparable intensities are used rather than 
one. The two control laser lead to a quasi-stationary 
probe-field pattern, known as stationary light [ij, llTI ]. 
We show that the dark polaritons of stationary light 
[9] behave like massive particles in 3D and can undergo 
Bose-condensation at hight temperatures. 
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FIG. 1: (Color online) a) The Raman interaction of two coun- 
terpropagating control lasers of Rabi- frequencies Q± and op- 
posite circular polarization coupling to the \s) — \e±) transi- 
tions of a double A system generates a quasi- stationary pat- 
tern of Stokes fields E±^ called stationary light. The station- 
ary dark-state polaritons formed in this interaction do not 
decay radiatively and have a quadratic dispersion profile in 
all directions, b) Off resonant coupling of the Stokes fields 
E± with other excited states leads to ac-Stark shift of lower 
level |s) resulting in a resonantly enhanced Kerr- type self- 
interaction of E±. 

We here consider a generalization of the original 
stationary-light scheme |14], which is shown in FigjT^). 
It involves two parallel A transitions consisting of the 
common ground levels \g) and \s) and the excited states 
|e±), coupled to opposite circular polarizations of quan- 
tized probe fields E± and control lasers described by the 
Rabi- frequencies Q.± [9|. Both A schemes are in two- 
photon resonance which guarantees EIT. For the present 



discussion we will assume that the control fields are ho- 
mogeneous and constant in time and thus set rt± = Q"^. 
The advantage of this set-up as compared to the original 
one, where both counterpropagating control laser cou- 
ple to the same dipole transition is that in the present 
scheme higher spatial harmonics of the ground-state co- 
herence do not occur and that there is no point in space 
where the total intensity of the control laser vanishes. 
The single-photon detunings of the upper states are de- 
noted as A+ and A_ respectively. The fields £^+ (^-) 
and Q-^ (^-) propagate in the -\-z {—z) direction with 
wavenumbers kp {—kp)^ and kc {—kc). 

We introduce normalized field amplitudes £± that 
vary slowly in space and time via E±(v^t) = 

'f- (^±(r, t) exp{— z(cjpt =f kpz)) + h.a.L and continu- 
ous atomic-fiip operators a^u{T,t) = ^ EjGAy(r) ^^^^ 
with a^^ = \/J^)jj{i^\ being the fiip operator of the jth 
atom. The sum is taken over a small volume AV around 
r containing AN atoms. 

Assuming that in the absence of the probe fields the 
control laser pump all population in the ground state |^), 
the dynamical equations read in the linear response limit, 
i.e. for a small probe intensity 
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Here n is the atom density and g = ^\/^ is the com- 
mon coupling constant of both probe fields with p denot- 
ing the dipole matrix element. T± = 7 + zA±, where 7 is 
the transverse decay rate of the transitions \e±) — \g) and 
kp = cOp/c is the carrier wavenumber of the probe field. 
If the characteristic time- and length-scales of the probe 
fields T and L are sufficiently large, such that the adi- 
abaticity conditions T <C {Labs/c, 7~^} and L ;» Labs 
hold, where Labs = ic/{g'^n) defines the resonant absorp- 
tion length, the system can be described by polariton-like 
quasi-particles [9(. One of these solutions, called dark- 
state polariton (DSP) does not involve excited states and 
is thus immune to spontaneous emission. 



^(r,t) = cos6> ('cos(/)^+(r,t) + sin(/)^_(r, t)]- 



-sin^S'(r, t) 

(4) 

where S = \fnc)gs (r, t) and the mixing angles are de- 
fined by tan2 6> = g'^n/Vfi with Q? = n\ ^ Q?_, and 
tan^ (j) = ^^_/^\. The longitudinal dispersion relation of 
the stationary DSPs is plotted in Fig |2] for tan(^ = 1 (red 
line). Also shown are the energies of the other eigensolu- 
tions. It proves useful for calculations to introduce spe- 
cial superpositions of these eigensolution termed bright- 



state polaritons. 

I>i(r,t) = -sin0£:+(r,t) + cos(/)^_(r,t), (5) 

l>2(r,t) = - sin 6>( cos (/)£:+ (r, t) + sin (/)£:_ (r,t)) 

-cosO S{Y,t). (6) 
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FIG. 2: (Color online) Dispersion of stationary dark-state 
polariton (red) and other quasi-particle solutions (blue) in 
propagation direction for equal detunings and control field 
Rabi- frequencies A+ = A_ = A, and Q+ = Q- = Q, 
tan^ = tan(/) = 1, A/Q = 2. Since only the real parts of 
the eigenvalues are of interest here, we set 7 = 0. One clearly 
recognizes a quadratic profile for the dark polariton around 
k = 0. It should be noted that the polariton is defined by 
the slowly varying envelopes of field and matter variables and 
thus the wavenumber k refers to spatial modulations of the 
slowly varying amplitudes. 

It is easy to derive the equation of motion of the DSP 
up to second order including the transverse directions. 
One finds for A± = A ::^ 7 and (j) = 7r/4 
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where Vgr = c cos^ 6 is the group velocity of the polariton 
if only one of the two coupling fields is present. As op- 
posed to microcavity exciton-polaritons, stationary DSPs 
behave as free, massive Schrodinger particles in all three 
spatial directions with a tensorial mass, m± = mv^Qc/'^gvi 
m7^^ = m^ 2/cpLabs (^/t + 0- Here we have introduced 
the recoil velocity of the atom at the probe frequency 
'^rec = hkp/m and the mass of the atoms m as compar- 
ative scales. The effective mass can be varied via the 
strength of the control laser, i.e. via the mixing angle ^, 
and is typically several orders of magnitude smaller than 
the mass of the atoms. The longitudinal mass my is again 
smaller than the transverse mass by the ratio of resonant 
wavelength to absorption length and the inverse normal- 
ized detuning 7/ A. It also has a small imaginary com- 
ponent, which describes the absorption of high-/c com- 
ponents well known from FIT [l6[. As can be seen from 



Fig|2]the quadratic dispersion of the dark- state polariton 
turns over into a linear dispersion at higher frequencies 
and eventuahy saturates giving rise to a finite band of fre- 
quencies whose width is of the order of the single-photon 
detuning A for |A| ^ ^,7- 

Using the well known expression for the critical tem- 
perature of condensation for a homogeneous, ideal gas 
and expressing the polariton mass in terms of the mass 
of the atoms one finds 



natom fPDSP_y^^V^ 
J Vr, 
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Here we have introduced the critical temperature of the 
atoms ksT^^'"'^ = 27rn'^/^ /[C{3/2)'^/^m] as a comparative 
scale, pdsp is the density of dark-state polaritons, which 
has to be smaller than the atomic density n in order to 
stay within the linear response limit. Since the group 
velocity Vgj- can be orders of magnitude larger than the 
recoil velocity of atoms Vj-ec and since kpL^^js ^ 1, the 
critical temperature can be much larger than that of the 
atoms. E.g. for Pdsp/^ = lO""*^, Vg^ = 1 km/s, v^qc = 5 
cm/s, Labs = 1 cm, A/7 = 10 and kp = 27r/500 nm, 
yields: Tc/T^^^^ ?^ 6 x 10^, i.e. a value in the mK regime. 
On the other hand the finite frequency window of EIT 
at large detuning ksT^iT ~ ^A limits the maximum al- 
lowed temperature where losses are small to the regime 
of 0.1 ... 1 mK. The losses outside the EIT transparency 
window may be used for evaporative cooling of polari- 
tons, which is however outside the scope of the present 
paper. 

Condensation is achieved differently from atoms or 
exciton-polaritons. Instead of changing the tempera- 
ture or the density one can dynamically change the criti- 
cal temperature of condensation by varying the effective 
mass, or respectively 6. Initial preparation of low-energy 
DSPs can be realized in various ways. Spontaneous Ra- 
man scattering cannot be used since it results in too high 
polariton momenta on the order of the recoil momentum 
corresponding to polariton temperatures on the order of 

Trec^^/ [m/^ mJ ). Direct RF coupling or storage of a 
coherent light pulse yield a low energy coherent polari- 
ton distribution. Alternatively storage of an incoherent 
light pulse results in the preparation of a non-equilibrium 
mixed state. 

Since the maximum temperature allowed by the fi- 
nite transparency window is below the temperature range 
where collisions in a gas cell could provide sufficiently fast 
thermalization, the latter should occur by elastic colli- 
sions between dark-state polaritons. These can be in- 
duced by resonantly enhanced optical nonlinearities pro- 
vided by the stationary-light coupling scheme itself. If 
e.g. the excited states in FigH] are hyperfine states then 
other off-resonant couplings exist, such as those shown in 
FigHb- This coupling induces ac-Stark shifts which give 
rise to a resonantly enhanced optical Kerr nonlinearity 



|l8l . Il9j . If we ignore off-resonant couplings of the con- 
trol laser fields, which is justified if the energy splitting 
between the lower states \g) and \s) is larger than |1^±|, 
the resulting effective Kerr-interaction reads 
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Axerr is the ouc-photou detuning from the additional ex- 
cited states, and " : :" denotes normal ordering. Since ([9|) 
only effects the electromagnetic component of the polari- 
tons, the optical Kerr interaction involves also nonlinear 
scattering of dark- into bright-state polaritons. Assum- 
ing sin^ ^ 1, i.e. Vgr <C c, one finds 
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^t|>2 ) ^. (10) 



The bright polaritons decay rapidly and thus the non- 
linear coupling of polariton modes results in an effective 
loss. Adiabatic elimination of the fast decaying bright 
polaritons yields in Born-Markov approximation a Liou- 
ville equation for the density matrix of the DSPs with 
a Hamiltonian part describing elastic two-body collisions 
and a nonlinear loss part. 
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From these one can extract the rate of elastic collisions 

Fcoll = 
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as well as the rate of collision induced losses 
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One notices that the characteristic time scale r of both 
processes is determined by the ratio of absorption length 
to group velocity. For typical vapor values such as Vgr = 1 
km/s and Labs = 1 cm one finds r = 10~^ s, while for 
a solid or a condensate of atoms where Labs can be as 
small as 10 /im one finds r = 10~^s. Since 7/AKerr 
as well as Pdsp/^ are small compared to unity, elastic 
collisions happen always on a much shorter time scale 
than the nonlinear losses. 

In addition to the nonlinear loss there is also a linear 
absorption resulting from the imaginary part of the po- 
lariton mass. It corresponds to the absorption of high 



spatial- frequency components in EIT. If L denotes the 
characteristic longitudinal length scale of the polariton, 
the linear loss rate can be estimated from eq.([7j): 
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(14) 



In order for the elastic collisions to be fast compared to 
the linear losses it is necessary that the optical depth 
OD = L/Ls,\)s of the medium over the length L of the 
polariton wavepacket fulfills 
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(15) 



Taking the above example of Vgr = 1 km/s and Labs = 1 
cm or 10 /im respectively, and assuming Axerr = IOO7 
and Pdsp/^ = 0.1 finally yields an elastic collision rate 
of Pel ~ 10^ s~^, or 10^ s~^. Thus for sufficiently large 
elastic collision rates one needs a large optical depth 
0D> 30, which is feasible however in solid-state systems 
or magneto-optical traps. 

The signature of condensation is here the macroscopic 
occupation of the k = momentum state. The tran- 
sition into modes with vanishing transverse momentum 
k^ = can easily be observed by switching off one of the 
two control laser. In this case the stationary polariton 
will be transformed into a moving one and will leave the 
sample as a light pulse. ^From the transverse profile of 
the emitted pulse one can deduce the population of the 
k^ momentum modes as indicated in Fig. [3l 
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FIG. 3: (Color online) toy: detection of condensation: switch- 
ing off one of the two control beams transforms the stationary 
polariton into a propagating polariton which can be observed 
as light pulse emitted from the sample, bottom: Transverse 
emission profile above (left) and below (right) condensation. 
An ideal quasi-homogeneous Bose-gas was assumed with fi- 
nite transversal extend modelled by a Gaussian distribution 
of width 178 Xt,x with Xt,x being the transversal de-Broglie 
wavelength at temperature T. 

In summary we have shown that stationary light po- 
laritons behave as massive Schrodinger-like particles with 



small and variable effective mass. In the presence of a 
resonantly enhanced optical Kerr nonlinearity the sta- 
tionary DSPs undergo elastic collisions at a rate deter- 
mined by the group velocity divided by the absorption 
length multiplied by the inverse normalized Kerr detun- 
ing. For a sufficiently high optical density of the sample 
the elastic collisions are fast enough to mediate a stim- 
ulated transition into the k = polariton mode. The 
critical temperature for this can be orders of magnitude 
larger than that of atoms. Condensation can be observed 
by releasing the stationary field into a propagating field 
and observing the transverse mode profile as well as the 
temporal coherence. In contrast to exciton-polaritons no 
cavity is needed and the polariton gas is three dimen- 
sional. As a consequence the stimulated transition into 
the k = mode is a true Bose-Einstein condensation. 



The authors would like to thank E. S. Demler, G. Mo- 
rigi and M. Weitz for stimulating discussions. The finan- 
cial support of the DFG through the GRK 792 and of 
the center of excellence OPTIMAS is gratefully acknowl- 
edged. 



[1] 
[2] 

[4] 
[5] 

[6: 

[7] 
[8] 
[9 

[10: 

[11 

[12 
[13 

[14 

[15 
[16 

[17 

[is: 

[19 



S.N. Bose, Z. Phys. 26, 178 (1924). 

A. Einstein, Sitz. Ber. Preuss. Akad. Wiss. (Berlin) 22, 
261 (1924). 

J. M. Blatt, K. W. Boer, and W. Brandt, Phys. Rev. 
126, 1691 (1962). 

P. Littlewood, Science 316, 989 (2007). 
J. Kasprzak, et a/., Nature 443, 409 (2006). 
R. Balili, V. Hartwell, D. Snoke, L. Pfeiffer, and K. West, 
Science 316, 1007 (2007). 

S.O. Demokritov et a/.. Nature 443, 430 (2006). 
D. Snoke, Nature 443, 403 (2006). 

F.E. Zimmer, J. Otterbach, R.G. Unanyan, B.W. Shore, 
and M. Fleischhauer, Phys. Rev. A 77, 063823 (2008). 
M. Fleischhauer, and M.D. Lukin, Phys. Rev. Lett. 84, 
5094 (2000). 

L.V. Hau, S.E. Harris, Z. Button, and C.H. Behroozi, 
Nature 397, 594 (1999) 

D.F. Phillips, A. Fleischhauer, A. Mair, R.L. Walsworth, 
and M.D. Lukin, Phys. Rev. Lett. 86, 783 (2001). 
C. Liu, Z. Button, C.H. Behroozi, and L.V. Hau, Nature 
409, 490 (2001). 

M. Bajcsy, A.S. Zibrov, and M.D. Lukin, Nature 426, 
638 (2003). 

S.E. Harris, Physics Today 50, 36 (1997). 
M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Rev. 
Mod. Phys. 77, 663 (2005). 

F.E. Zimmer, A. Andre, M.D. Lukin, and M. Fleis- 
chhauer, Opt. Comm. 264, 441 (2006). 
H. Schmidt, and A. Imamoglu, Opt. Lett. 21, 1936 
(1996). 

D.E. Chang et oL, Preprint at 

( |http://arxiv.org/abs/0712.1817| ) (2007). 



